Detection and engineering of spatial mode entanglement with ultra-cold bosons 
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We outline an interferometric scheme for the detection of bi-mode and multi-mode spatial entan- 
glement of finite-temperature, interacting Bose gases. Whether entanglement is present in the gas 
depends on the existence of the single-particle reduced density matrix between different regions of 
space. We apply the scheme to the problem of a harmonically trapped repulsive boson pair and 
show that while entanglement is rapidly decreasing with temperature, a significant amount remains 
for all interaction strengths at zero temperature. Thus, by tuning the interaction parameter, the 
distribution of entanglement between many spatial modes can be modified. 
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I. INTRODUCTION 

Understanding and controlling entanglement in many- 
body systems is one of the most important challenges in 
quantum mechanics today. The rewards are significant 
and are expected to not only lead to new insights into the 
properties of solid state systems and phase transitions 
but also to new designs for highly efficient quantum 
information devices. Ultra-cold bosonic gases offer an 
ideal arena to explore many-body entanglement, as ex- 
perimentalists have at their disposal designer condensed 
matter systems whose parameters can be controlled with 
unprecedented precision . 

Entanglement often exists naturally in the ground 
state of a many-body system Q , where it resides between 
the degrees of freedom of the particles and is a property of 
the first quantised many-body wavefunction. However, in 
ultra-cold gases the particles are inherently indistinguish- 
able, which requires the symmetization of their many- 
body wavefunction and means that the Hilbcrt space no 
longer has the tensor product structure required to define 
entanglement. The first quantised many-body wavefunc- 
tion of indistinguishable particles may therefore contain 
quantum correlations but such correlations are usu- 
ally considered unable to violate a Bell inequality or pro- 
cess quantum information [3, S 0] • 

Ultra-cold gases are also well described within the 
framework of second quantisation, where instead of work- 
ing directly with the many-body wavefunction, one de- 
fines a complete set of field modes that are occupied by 
particles. Second quantisation therefore offers the possi- 
bility of entanglement between modes. Entanglement is 
dependent on the choice of modes, but provided the cor- 
rect choice is made, investigating entanglement between 
distinguishable modes 0, H circumvents the difficul- 
ties of defining entanglement between indistinguishable 
particles 
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To illustrate the differences between particle and mode 
entanglement, let us consider two non-interacting bosons 
in a trap at zero temperature. In first quantisation, the 
wavefunction is the symmetrized product, ^'12 (a;, y) = 
^(0i(a;)02(y) 4- (l)i{y)(f>2ix)), where (^(x) is the ground 
state of the confining potential. No entanglement ex- 
ists between the particles, since indistinguishability for- 
bids us from assigning to any particle a specific set 
of degrees of freedom. Conversely, in second quanti- 
sation one can define a pair of spatial modes, A and 
B, where each mode occupies half the confining geom- 
etry. Since both the particles are coherently distributed 
over these modes, the system is described by the en- 
tangled state, li^As) = 5(120) + V2\ll) + |02)), where 
\mn) = \m)A ® denotes m particles in mode A and 
n particles in mode B (with m + n ~ 2). 

In this paper we outline a scheme for the detection of 
bi-mode and multi-mode spatial entanglement for a fi- 
nite temperature, interacting Bose gas of any (including 
unknown) particle number. We show that entanglement 
is detected via the single-particle reduced density matrix 
(SPRDM). We apply our scheme to the example of a 
harmonically trapped, interacting boson pair [l2|, where 
the SPRDM also acts as a quantifier of entanglement. 
We find that for all interaction strengths, entanglement 
between pairs of modes rapidly decreases with tempera- 
ture. While at zero temperature a significant amount of 
entanglement remains even in the limit of infinite inter- 
action. Moreover, we note that our detection scheme is 
also relevant to recent proposals to observe non-locality 
of single particle between spatial modes [l^, [1^1 ■ 



II. DETECTION SCHEME 

The correlations of entanglement are locally basis in- 
dependent, so that one needs to measure each mode in at 
least two bases in order to differentiate them from classi- 
cal correlations. While a superselection rule that forbids 
coherent superpositions of cigenstatcs of different mass 
seems to rule such measurements out for any atomic 
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FIG. 1: Schematic showing a trapped wavefunction spht into 
two spatial modes, a and b. The modes are combined at a 
50/50 beamsphtter and the particles in the output modes c 
and d are counted. 



system, recent work has shown that such measurements 
are theoretically possible Moreover, it has recently 
been predicted that this natural mode entanglement of 
massive particles can be used as a resource for quantum 
communication [l6j . However, in practice it will be diffi- 
cult to locally manipulate, i.e. rotate, the spatial modes. 
In the following we will show that spatial entanglement 
can also be detected and quantified by making global op- 
erations on the modes. 

Let us consider a gas in a confining geometry (sec 
Fig.[T]) which is mathematically, but not necessarily phys- 
ically, divided into two, non-overlapping spatial modes, 
a and b. The field operators, i^l = j^dx g*{x)'ii)^ {x) and 
ipi = dx g(x)ip{x), create and destroy particles in mode 
i = a, b, where \g{x)\'^dx = 1 ensures that the commu- 
tation relations, [■(/',;, V^j] = Sij, are satisfied. The quan- 
tity, g{x), specifies how the set of points in a spatial mode 
are averaged over. To generate interference, the parti- 
cles in the two spatial modes, a and b, are mixed at a 
50:50 beamsplitter, which transforms the input modes as 
^1 = 7l(ViJ + ^t) and = ^ii'l - After the 

beamsplitting operation, the number of particles in the 
output modes, c and d, are counted and compared to the 
fully separable case. If the number of coincidences is dif- 
ferent to the separable case, we can conclude that there 
must have been entanglement between the spatial modes. 

A bosonic gas of fixed particle number, N which is in 
a fully separable state w.r.t. the spatial modes a and b 
can be written as 



N 



ra=0 



where ^„ p„ = 1. The beamsplitter transforms this state 
such that if the total number of particles is even, one 
detects the same number of particles in each of the output 
modes, i.e. A7V = \Nc — = 0. For an odd number 
of particles an ensemble average leads to the same result 

Conversely, if the initial state of fixed particle number 
is of an arbitrary form, p, the difference in particle num- 
bers detected in the modes c and d may be non-zero due 
to entanglement between the modes. 
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(2) 



where 



eab = tr 



dx / dx'g{x)g*{x')p^'^\x,x'). 



Pn\n){n\a(E>\N -n){N -n\b, (1) 



(3) 

We have derived the above result using the Fourier de- 
composion of the field operators, ip^ {x) and ipi^) in terms 
of the momentum modes, (f>k{x), as ip\x) = 4>X{x) 
and likewise for tp{x). When tab is non-zero, the state p is 
different to the separable case and is therefore necessarily 
entangled w.r.t. the bi-modal split into a and b. 

The quantity Eab is given by the off-diagonal ele- 
ments of the SPRDM which is defined as p^^'>{x,x') = 
ni:4)j:{x)(f>*^{x'). Here rij: is the number of bosons that 
occupy the k-th momentum mode. The SPRDM is a one- 
body correlation function, which characterises important 
coherence properties of a many-body system [l8| and its 
off-diagonal elements are related to the visibility of inter- 
ference fringes in a two slit experiment . 

Our scheme therefore applies to systems whose correla- 
tions from a basic group of one particle [lO] ■ States of the 
form, 1 20) + |02), whose correlations are second order are 
not detected by the SPRDM. Such states require careful 
engineering of bosonic systems [2l|, although they may 
form naturally in fcrmionic systems, i.e. Cooper pairs in 
superconductors. Our scheme is applicable to all bosonic 
gases trapped in orthodox geometries. 

Next, let us show that eab can also be used to quantify 
mode entanglement for certain systems, one of which is 
the boson pair model discussed below. For this eab must 
fulfil three basic criteria (2^ : (i) is zero when the 
state is separable, (ii) eab is invariant under local uni- 
tary operations and (iii) eab docs not increase under lo- 
cal general measurements and classical communication 
(LGM+CC). The validity of (i) is shown above and (ii) 
is guaranteed since the trace is basis independent. One 
can prove (iii) as follows. To implement LGM-I-CC, the 
two spatial modes, a and 6, are each coupled to a local 
environment by a general completely positive map. The 
environments are allowed to communicate classically ad 
infinitum and the total number of particles in the gas and 
environment is fixed. One can then show that e^h does 
not increase under LGM-I-CC. 

The above scheme can also analyse multi-mode en- 
tanglement, i.e. the simultaneous entanglement of more 
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than two modes. In the following we use the differ- 
ent notions of separability discussed in [2^. A general 
M-mode state is fully separable and contains no entan- 
glement if it is a convex combination of states for each 

mode, Psep(M) = J2iPiPi^^ (8) • • • (8) p-*^^ Here each com- 
posite state, p~ , corresponds to a single spatial mode 
with a fixed number of particles. On the other hand, 
entanglement may be present between certain subsets of 
spatial modes, for which one can define a p-separable 

state, Psepip) = Y.iPiP'P ® ■■■ ® p'f\ where p < M 
with equality when the state is fully separable as above. 

When the composite state, , describes more than one 

spatial mode, the spatial modes contained within pp' 

are necessarily entangled otherwise pf would be writ- 
ten as a product of states for the individual modes, 
= ppi) (g) • • • (g) /Sp"' «)•••, and the overall state of the 
system would be cr-separablc, where p < a < M [2^ . 

To determine whether a given M-mode state contains 
multi-mode entanglement, one can check a bi-partite en- 
tanglement criterion between all 2^^^^ — 1 unique divi- 
sions of the system into two blocks, A and B. Depending 
on which pairs of blocks are found to be separable one 
can conclude that the state has entanglement between 
different subsets of spatial modes. Here the bipartite 
entanglement criterion is eab = tr[^^^B/5] of eq. 
i.e. the SPRDM between two blocks of spatial modes, A 
and B. The field operators, ^'^ and ^x, for blocks of 



spatial modes are defined as VP^j- = ^ 
J2iex where X = A,B and Y,i 

that the commutation relations, 
satisfied. 
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Ci -ipl and *A- = 
c, P ~ 1 ensures 



= Sx^Y, are 



We now describe how entanglement changes for the (i) 
fully separable (ii) the p-separable and (iii) the fully en- 



tangled states: (i) The fully separable state, p 



sep{M) 



ad- 



mits no interference between all 2^^~^ — 1 partitions into 
blocks, i.e. eAB = for all A and B. (ii) A p-separable 
state, Psep{p) where p < M, admits no interference for 
2P~^ — 1 partitions into blocks from the total of 2^^~^ — 1 
partitions, i.e. eAB ~ for 2^^^ — 1 choices of A and B. 
(iii) A fully entangled state, p = ^^PiPi, admits inter- 
ference for all 2^^"^ — 1 choices of blocks, i.e. eab 
for all A and B. The fully entangled state, p, necessarily 
contains some form of multi-mode entanglement, since 
otherwise there would exist a decomposition such that p 
is p-separable. 



III. BOSON PAIR MODEL 

In the following we will apply our scheme to the phys- 
ically realistic model of a harmonically trapped pair of 
ultracold, interacting, bosonic atoms in effectively one 
dimension. The Hamiltonian of such a system is given 
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FIG. 2: Amount of entanglement between two spatial regions 
occupied by a boson pair as specified in the text as a function 
of temperature. The values of the dimensionless interaction 
parameter increase from the top curve to the bottom curve as 
giD = 0, 2, 5, 10, cxD. Temperature is scaled in units of hio/ks- 
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H 



E 



2 dxj 



+ 7:x^,] +9idS{\x,-Xj\), (4) 



where all lengths are scaled in units of the ground 
state size and all energies in units of the harmonic fre- 
quency. The one dimensional coupling constant, gio, 
is related to the three dimensional s-wave scattering 
= (a± - Ca3Dr\ where C is the 

maj_ ^ oj-y / 1 

1.4603 UM- This Hamiltonian can be 



length by gm 
constant, C — 



decoupled by moving into the centre of mass and rel- 
ative coordinate frames labelled by X and x, respec- 
tively, [T^ and the two-body wavefunction can subse- 
quently be written as n,v{x\,X2) = '4'n{X)ip^(x). The 
eigenvalues for the centre of mass motion are given by 

= 0, 1 ... , with corresponding 

MnHn{X)e ~. Here Mn is 
the normalisation constant and Hn{X) are the Hermite 
polynomials. For the relative motion, the single parti- 
cle eigenstatcs are '4}]f^[x) = Mv e~^U (|; — 
where v — 0,2,4... and the U{a,b,z) are the con- 
fluent hypergeometric functions. The corresponding 
eigenenergies, E,^, are determined by the roots of the 



= (n -f i) for n 
eigenstates, ip^°"'{X) = 



implicit relation —giD 



h2 



The 



'r(-^+i)^ 

kernel of the density operator in position representa- 
tion is Pnv{x,x' ,X2) — ^n_^{x,X2)^n,y{x' ,X2)j with the 

SPRDM defined as pni>{x,x') = Pnuix, x' , X2)dx2 ■ 
The SPRDM in thermal equilibrium is given by 
p'^^\x,x') = J2'i^J27PnuPnuix,x'), whcrc Pnu is the 
Boltzmann weight, Pm, = exp(^^^r-), k), is the Boltz- 
mann constant and Z is the partition function. 
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FIG. 3: The system is split into three spatial modes, a,b and 
c. Mode a is defined as the region — oo < x < —i^ i mode b as 
— < a:: < and mode c as < x < cx). Entanglement is 
investigated between neighbouring modes, a and 6 (b and c) 
(left hand side) and also between the outer two modes, a and 
c (right hand side) at T = as interaction is varied. Three 
different central mode lengths, L^, = 4, 2.8, 1.4 are taken. 



IV. BI- AND MULTI-MODE ENTANGLEMENT 
RESULTS 

We define two equal length modes, a and 6, and calcu- 
late their spatial entanglement for the boson pair model 
as a function of temperature and interaction strength us- 
ing cq. ([3]). The results are displayed in Fig. [21 With 
increasing temperature the gas becomes a statistical mix- 
ture of momentum modes, which destroys the 
fixed relative phase between the spatial modes and con- 
sequently the entanglement between them is severely de- 
pleted. The presence of interaction also degrades the 
quality of the entanglement at T = 0, but a significant 
amount persists even in the Tonks- Girardeau limit, of 
impenetrable bosons, qid = oo. At first this may seem 
surprising, since in the Tonks- Girardeau limit the system 
can be mapped onto an ideal fermionic atom pair [2^, 
for which one may not expect any entanglement to be 
present. However, only the local properties of the bosons 
become identical to free spin-polarised fermions due to 
the Bosc-Fcrmi mapping and therefore entanglement is 
not affected. 

We also examine pairwise entanglement in this model. 
For this space is divided into three modes, a, h and c, 
and we check entanglement between all pairs, ah, he and 
ac, as a function of the interaction strength at zero tem- 
perature. The results are shown in Fig. [31 where three 
different central mode lengths, L;,, are considered. The 
left graph of Fig. [3l shows entanglement between neigh- 
bouring modes, ah (and for symmetry reasons he). Here 
entanglement decreases as the interaction strength, .giu, 
increases for all lengths studied. The right hand side 
graph of Fig. [3l shows the entanglement between outer 
modes, a and c. The entanglement increases initially for 
all mode sizes as the interaction strength is varied. For 
the central mode size of if, = 1.4 entanglement how- 
ever is not a monotonic function of giu and decreases 
after the initial rise. For Lf, = 2.8 and 4 this is not the 
case. The behaviour of entanglement in both graphs can 
be explained by the dclocalisation of the bosons away 



from the centre of the trap due to greater interaction 
strengths. Stronger interaction therefore increases the 
correlation length (i.e. the entanglement between the 
outer two modes), yet decreases the overall coherence of 
the sample. 

We also check the existence of multi-mode entangle- 
ment in our model. For this we have divided our two- 
particle wavefunction into M modes and calculated ^ab 
for the 2^^^^ — 1 unique partitions of these modes into 
two blocks, A and B. Recall that, eab, must be non- 
zero for all possible block combinations for the state to 
be fully entangled. We have found that at zero tempera- 
ture the system possesses full multi-mode entanglement, 
i.e. eAB 7^ for all A and B, for a wide range of values of 
the interaction parameter, provided that the set of modes 
are defined within the coherence length of the sample. 

Here we have used the repulsively, interacting boson 
pair model to illustrate our scheme, since it is analytically 
tractable and contains all the essential features of larger 
models. Our scheme to detect bi- and multi-mode entan- 
glement can also be applied to Bosc gases with a greater 
number of particles (for all interaction strengths and tem- 
peratures), but qualitatively the results will remain un- 
changed while the task of computing the single-particle 
reduced density matrix between modes will quickly be- 
come very demanding. We note that bi-mode entangle- 
ment of a non-interacting Bose gas has been studied be- 
fore for N particles at zero temperature and at finite 
temperatures [26| . 

Our scheme can be implemented using currently avail- 
able technologies. Atomic beamsplitters can be realised 
with optical potentials [27[ and separate modes can be 
defined using spatially selective outcoupling techniques 
[l^. To perform beamsplitting, on systems with strong 
interactions one should raise a potential barrier between 
the desired modes on a non-adiabatic time scale so that 
the coherences between the regions are not lost, (i.e. so 
that the gas does not enter a Mott-like state). The poten- 
tial barrier should then be lowered rapidly enough so that 
the coupling between the wells is greater than the on-site 
interaction energy of the wells. For instance, one could 
envisage a setup where the trap is swiftly modulated from 
harmonic to double-well potential and is then switched 
off and the sample left to interfere. The entanglement 
can be inferred according to eq. ^ from the visibility of 
the resulting interference fringes, see for example p^.[28j. 
where mode entanglement between regions of space has 
already been measured indirectly. 



V. CONCLUSIONS 

We have outlined a scheme that detects bi-modal and 
multi-modal spatial entanglement of cold bosonic gases 
using simple atom-optic techniques. We show that spa- 
tial entanglement is detected by the SPRDM between 
different modes, a quantity that is related to the visibil- 
ity of routinely measured interference fringes. We have 
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demonstrated our scheme using the model of a harmoni- 
cally trapped boson pair and have found the existence of 
bi- and multi-mode entanglement within the gas. For all 
interaction strengths, increasing temperature rapidly de- 
grades the amount of entanglement, but at zero temper- 
ature entanglement still remains, even in the presence of 
strong interactions. Therefore, we have shown that the 
ability to vary the particle interaction strength allows 
one to engineer the distribution of entanglement over the 



trap, thus allowing a tunable source of entanglement. 
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